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Chapter 1. SAMPLE WAVE PROBLEMS 

To describe a problem in mathematical terms, one must make use of the basic laws 
that govern the elements of the problem. In continuum mechanics, these are the conser- 
vation laws for mass and momentum. In addition, empirical constitutive laws are often 
needed to relate certain unknown variables; examples are equations of state, Hooke's 
law between stress and strain, etc. 

To derive the conservation law one may consider an infinitesimal element (a line 
segment, area or volume element), yielding a differential equation directly. Alternately, 
one may consider a control volume (or area, or line segment) of arbitrary size in the 
medium of interest. The law is first obtained in integral form; a differential equation is 
then derived by using the arbitrariness of the control volume. The two approaches are 
completely equivalent. 

Let us first demonstrate the differential approach. 



1 Transverse vibration of a taut string 

Referring to Figure 1, consider a taut string stretched between two fixed points at a; = 
and X — L. Let the cross-sectional area be S. If there is an initial stretching of AL, the 
initial tension T must be 

AL 

T^ES— 
ij 

by Hooke's law, where E is Young's modulus. 

Now study the lateral displacement of the string from the initial position. By the law 
of conservation of transverse momentum, the total lateral force on the string clement 
must be balanced by its inertia. Let the lateral displacement be V[x, t) and consider a 
differential element between x and x-\-dx. The net transverse force due to the difference 
of tension at both ends of the element is 

(Tsina)^^^^ - (Tsina)^, 
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Figure 1: Deformation of a taut string 



where 



dV 



sma — 



dV 
dx 



We shall assume the displacement to be small everywhere so that the slope is also small: 
1^ <S 1. The local value of sin a can then be approximated by 

dV ^ fdV 

dx \ dx 

where the expression 0{5) stands for of the order of S. For any smooth function /, 
Taylor expansion gives 



f{x + dx)-f (x) 



-J- ] dx + O idx)^ , 
ox ' 



where the derivative is evaluated at x. Hence the net tension is 
The instantaneous length i{x, t) of the string from to x is 



i{x,t) 



dx 



^ ~^ ^ dx^ 



1/2 



^ ^ ^ ^dx 



It follows that 



^ ""^of^^ for all 0<x<L, 



X 



which is of second-order smallness. The string length, hence the tension, is essentially 
unchanged with an error oiO{dV/dxy, i.e., T can be taken as constant with a similarly 
small error. Thus the net tension in the string element is well represented by 

T—rdx. 
ox^ 
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If the mass per unit length of the string is p, the inertia of the element is p{d'^V/dt^)dx. 
Let the applied load per unit length be p(a;, t). Momentum conservation requires that 

pdx^r— = T—-^dx + pdx + O idx) . 
ot^ ax^ 

Eliminating dx and taking the limit of dx — > 0, we get 

p d'^V d'^V p 



T df^ dx'^ T' 



(1.1) 



This equation, called the wave equation, is a partial differential equation of the second 
order. It is linear in the unknown V and inhomogeneous because of the forcing term on 
the right-hand side. 

Is the longitudinal displacement U important in this problem? Conservation of 
momentum in the x direction requires that 

P^^-Q^ = (Tcosa)^^^ - {Tcosa)^. 

Since 

dx 1 ^ fdV 

cos a = = , ^ 1 + O — 

the acceleration is of second-order smallness 

pd'^U ^( d {dV\^\ ^ f fdV\ p d^V 



T dt^ \dx\dx J \\dx J T df^ 



Hence U = O (^^^V hy twice integration with respect to t, and the longitudinal dis- 
placement can be ignored. 

The differential equation (1.1) involves second-order derivatives with respect to both 
X and t. Two auxilliary conditions are needed for each variable. For example, at the 
initial instant, we may prescribe both the displacement and the velocity: 

V{x,0)^f{x) (1.2) 

and 

dV 

— {x,0)^g{x). (1.3) 
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These statements are called the initial conditions. In addition we must also specify the 
boundary conditions at the ends. For a string stretched between two fixed ends, we 
require 

V{0,t)=0 and V{L,t) = 0. (1.4) 

Together with the partial differential equation, these auxilliary conditions define the 
initial-boundary-value problem. Prom the mathematical point of view, it is important 
to establish whether such a problem is well posed. This question involves the proof for 
the existence, uniqueness and stability of the solution. 

As seen in this example, Taylor expansion is used at almost every step of the deriva- 
tion. Indeed, it is indispensable not only in deriving governing equations, but also in 
obtaining approximate solutions of the equations, and in analyzing the physical content 
of the solution. 

Note that the dimension of the coefficient T/ p is 

ML/t^ fL^^ 



[velocity] 



2 



M/L \t 

Now introduce the notation c = y^T/p , which is a characteristic velocity of the physical 
problem. Equation (1.1) can then be written 

1 d^V d^V p 



c2 dt^ dx^ T' 
which is called the wave equation arising in numerous contexts. 



:i.5) 



2 Longitudinal vibration of an elastic rod 

Consider an elastic rod with the cross-sectional area S{x) and Young's modulus E, as 
shown in Figure (2). Let the longitudinal displacement from equilibrium be U{x,t). 
The strain at station x is 

AU dU 
lim — — = — — . 
Ax-^o Ax ox 

By Hooke's law, the tension at x is 
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U U+dU 








X x+dx 



Figure 2: Longitudinal deformation of an elastic rod 
Now the net tension on a rod element from x to x + dx is 



ES 



dx 



x+dx 



Let the externally applied longitudinal force be f{x,t) per unit length. Momentum 
conservation requires that 

pb-—T-dx = — — dx + tdx + (J{dx) . 

ot^ ox \ ox J 

In the limit of vanishing dx, we get the differential equation: 



(2.1) 



In the special case of uniform cross section, S = constant, and zero external force, 
then U satisfies the inhomogeneous wave equation 

1 d^U d^U f 



c2 dt^ dx^ ^ ES' 



(2.2) 



where c = \^Ejp has the dimension of velocity. 

The simplest boundary conditions are for fixed or free ends. If both ends are fixed, 
then, 

U{0,t)=0 and U{L,t) = 0. (2.3) 
If the left end is fixed but the right end is free, then 



U{0, t) = and 



dU 
dx 



{L,t) = 0, 



(2.4) 
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since the stress is proportional to the strain. Again, the most natural initial conditions 
are 

dU 

U{x,0) = f{x), and —{x,Q)=g{x), (2.5) 

where / and g are prescribed functions of x for < x < L. 
Let us change to the integral approach in the next example. 



3 Traffic flow on a freeway 

One of the mathematical models of traffic flow is the hydrodynamical theory of Lighthill 
and Whitham (1958). It is a simple theory capable of describing many real-life features of 
highway traffic with remarkable faithfulness. Consider any section of a straight freeway 
from X — ato X — b, Figure 3. Assume for simplicity that there are no exits or entrances, 
and all vehicles are on the go. Let the density of cars (number of cars per unit length of 
highway) at x and t be p{x,t), and the flux of cars (number of cars crossing the point 
X per unit time) be q{x,t). By requiring that the number of cars within an arbitrary 
section from a to fo be conserved, we have 

d- f'' 

~dij = ■ 

Rewriting the right-hand side 

cb 



q[b,t)-q[a,t)^ j g^dx , 



^' + ^)dx^O. (3.1) 



we get 

Ja ^dt ' dx 

Since the control interval (a, b) is arbitrary, the integrand must vanish. 

This result can be argued by contradiction, which is a typical reasoning needed to change 
an integral law to a differential law. Suppose that the integrand is positive somewhere 
within (a, 6), say, in the range (a', 6') G (a, 6), and zero elsewhere in (a, 6), then the 
integral in (1.3.1) must be positive. But this is a contradiction. The assumption that 
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Figure 3: (a). A section of the freeway, (b). The relation between traffic flux rate and 
traffic density. 

the integrand is positive somewhere is therefore wrong. By a similar argument, the 
integrand cannot be negative anywhere, and hence must be zero everywhere in (a, b). 

Equation (1.3.2) is the law of conservation of cars. Having two unknowns q and p, a 
constitutive relation between p and q is needed and must be found by field measurements. 
Heuristically, q must be zero when there is no car on the road, and zero again when the 
density attains a maximum (bumper-to- bumper traffic), hence the relation between q 
and p must be nonlinear 

q^q{p) (3.3) 
as sketched in Figure 1.4.b. With this relation, (1.3.2) becomes 

This result is a first-order nonhnear partial differential equation and will be used to 
deduce a variety of interesting phenomena of traffic flow. 

In all the examples studied so far the final governing equation involves only one 
unknown. Now we will examine a problem with several unknowns. 

4 Wave propagation in arteries 

We shall examine the pulsating flow of blood in an artery whose wall is thin and elastic. 
As a first exercise let us assume that there is only pulsation but no net flux. Because of 
the pressure gradient in the blood, the artery wall must deform. The elastic restoring 
force in the wall makes it possible for waves to propagate. 
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The artery radius a{x, t) varies from the constant mean Ug in time and along the 
artery (in x). Let the local cross sectional area^ he S — tto^, and the averaged velocity 
be u{x,t). Consider a fixed geometrical volume between x and x + dx, through which 
fluid moves in and out. Conservation of mass requires 



dS d{uS) 
dt dx 



0, 



(4.1) 



Next the momentum balance. The time rate of momentum change in the volume must 
be balanced by the net influx of momentum through the two ends and the pressure force 
acting on all sides. The rate of mometum change is 

d{puS) 



The net rate of momentum influx is 

djpu^S) 
dx 



dt 



, duS ^du 

ax — —pu— pub — 

dx dx 



(4.2) 



(4.3) 



The net pressure force at the two ends is 

d{pS) 



dx 



^dp dS 
dx ^ dx 



while that on the sloping wall is 



27rap 



da dS 
P- 



The sum of all pressure forces is 



dx dx 

,dp 
dx 



(4.4) 



Balancing the momentum by equating (4.2) to the sum of a(4.3) and (4.4) we get, after 
making use of mass conservation (4.1), 



/ du 



du\ 
dx J 



dp 
dx 



(4.5) 



Let the pressure outside the artery be constant, say zero. The change in the tube 
radius must be caused by the change in blood pressure. Refering to Figure 4, the elastic 
strain due to the lengthening of the circumference is 27rda/27ra = da/a. Let h be the 



^Denoted by A in Figure 4. 
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Figure 4: Forces on the artery wall. Note: The symbol for the cross- sectional 
area should be changed from A to S. 

artery wall thickness, assumed to be much smaller than a, and Young's modulus E. The 
change in elastic force is 2Ehda/a which must be balanced by the change in pressure 
force 2a dp , i.e., 

2Eh da 

= 2a dp, 



which implies 



dp Eh dp \fKEh , , 



Pressure increases with the tube radius, but the rate of increase is smaller for larger 
radius. Upon integration we get the equation of state 



p-Po = -E{h/a) = -x/ttE/i/v^. (4.7) 



Eq (4.5) may now be rewritten as 



where C is defined by 
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and has the dimension of velocity. In view of (4.6), equations (4.1) and (4.8) are a pair 
of nonhnear equations for the two unknowns u and S. 

Linecirization: 

For infinitesimal amplitudes we can linearize these equations. Let a — ao + a' with 
a' <^ Oo then the (4.1) becomes, to the leading order. 



The linearized momentum equation is 



The linearized form of (4.6) is 
which can be used in (4.11) to get 



du dp 



dp = ^da' (4.12) 

«0 



du Eh da' ■,^\ 



Finally (4.2) and (4.8) can be combined to give the wave equation: 



dt^ " dx^ 
where 

' Eh 



= cf-^ (4.14) 



2pao 

Alternately one can eliminate a to get an equation for u 



(4.15) 



d'^u ^d'^u /. -.ox 

Because of (4.12), the dynamic pressure is governed also by 

^ = (4 17) 

df' "dx^ ^ ^ 

All unknowns are governed by the same equation due to linearity and the fact that all 

coefficients are constants. 

Comments on linearization: 
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To find out the accuracy of linearization, it is useful to estimate first the scales of 
motion. Let A,T,L,U and P denote the scales of a',t,x,u and p' respectively. It is 
natural to take L = CoT. Prom (4.1), (4.6) and (4.5) we get the relations among the 
scales of dynamical quantities 

"^ = ^, hence U = 



EhA 



U IP 11 EhA 



It follows that 



T pL pL/T al 
A L 11 EhA 



hence. 



Oo p L al 

I^_Eh_^ 
T^~'p^o~''" 

With these scales the ratio of a typical nonlinear term to a linear term is 

4| U'/L^ U ^A 
f U/T L/T ao 

Hence the condition for linearization is that 

an 



5 Shallow water waves and linearization 

5.1 Nonlinear governing equations 

If water in a lake or along the sea coast is disturbed, waves can be created on the surface, 
due to the restoring force of gravity. Consider the basic laws governing the motion of long 
waves in shallow water of constant density and negligible viscosity. Referring to Figure 
5, let the z axis be directed vertically upward and the x, y plane lie in the initially calm 
water surface, h{x, y) denote the depth below the still sea level, and C,{x, y, t) the vertical 
displacement of the free surface. Take the differential approach again and consider the 
fluid flow through a vertical column with the base dxdy. 
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First, the law of mass conservation. The rate of volume increase in the column 



dt 



dxdy 



must be balanced by the net volume flux into the column from all four vertical sides. In 
shallow water, the horizontal length scale, characterized by the wavelength A, is much 
greater that the vertical length h. Water flows mainly in the horizontal planes with 
the velocity u(x, y, t), which is essentially constant in depth. Through the vertical sides 
normal to the x axis, the difference between influx through the left and outflux through 
the right is 

-[u{C + h) \^+ax - u{C + h)\^]dy ^ - ^^[u{C + h)] + 0{dx) | dxdy. 

Similarly, through the vertical sides normal to the y axis, the difference between influx 
through the front and outflux through the back is 

- (C + h) \y+dy - V {C, + h)\y]dx = - ^^[v {C, + h)] + 0{dy) | dydx. 

Omitting terms of higher order in dx, we invoke mass conservation to get 

-^dxdy = - I ^ (C + ^)] + ^ b (C + h)] + 0{dx, dy) j dxdy. 



In the limit of vanishing dx, dy, we have, in vector form. 



dt 



+ V • [u(C + h)] = 0. 



(5.1) 



This equation is nonlinear because of the quadratic product of the unknowns u and (. 



O 



z=-h 



■ / y / / / // /' 



dx 



dy 



Figure 5: A column element of fluid in a shallow sea 
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Now the law of conservation of momentum. In shallow water the vertical momentum 
balance is dominated by pressure gradient and gravity, which means that the distribution 
of pressure is hydrostatic: 

p = pg{C-z) , (5.2) 

where the atmospheric pressure on the free surface is ignored. Consider now momentum 
balance in the x direction. The net pressure force on two vertical sides normal to the x 
direction is 



dxdy 
ox 



/ pdz = -dxdy— / pg{C,-z)dz 
= -pg{C + h) ^^^^ ' dxdy. 



The hydrodynamic reaction from the sloping bottom to the fluid is 

-p—dxdy = pg{C + h)—dxdy. 

The change of fluid momentum consists of two parts. One part is due to the time rate 
of momentum change in the water column 

d 



[P^*(C + h)] I dxdy, 

and the other is due to the net flux of momentum through four vertical sides: 

d d 
^[P«^(C + h)]dxdy + —[puv{C + h)]dydx. 

Equating the total rate of momentum change to the net pressure force on the sides and 
on the bottom, we get 

I WC + h)] + ^[pu'iC + h)] + ^[puviC + h)] 

The left-hand side can be simplified to 

du du du\ ,^ ,, { dC d , d 

\{C + K) + ou<^ 
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by invoking continuity (5.1). Hence the x momentum equation reduces to 

^+u—+v—^-— (53) 
dt ^ dx ^ dy ^ dx' 

Similarly, momentum balance in the y direction requires 

dv dv dv d( 
dt ^ dx ^ dy ^ dy' 

These two equations can be summarized in the vector form: 

dw. 

— + u • Vu = -^VC- (5.5) 

Equations (5.1) and (5.5) are coupled nonlinear partial differential equations for three 
scalar unknowns u and C,. 

Now the boundary and initial conditions. On a shoreline there can be no normal 
flux, therefore, 

/iu-n = on S, (5.6) 

where n denotes the unit normal vector pointing horizontally into the shore. This 
condition is applicable not only along a cliff shore where h is flnite, but also on a 
shorehne where /i = 0, as long as the waves are gentle enough not to break. In the latter 
case the whereabout of the shoreline is unknown a priori and must be found as a part 
of the solution. 

At the initial instant, one may assume that the displacement C(a;,y,0) and the ver- 
tical velocity of the entire free surface ^C{x,y,0) is known. These conditions complete 
the formulation of the nonlinear shallow water wave problem. 



5.2 Linearization for small amplitude 

For small amplitude waves 

where A is the characteristic amplitude. Equation (5.1) may be simplified by neglecting 
the quadratic term 

+ V ■ hu = 0. (5.8) 

at 
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Denoting the time scale by the wave period T and the horizontal length scale by the 
wavelength A, we equate the order of magnitudes of the remaining two terms above to 
get 

A uh . ^ . A uT ^ 
- ~ — , implymg ^ 

Now let us estimate the importance of the quadratic term u • Vu in the momentum 
equation by assessing the ratio 



at 



A 



Clearly the quadratic term representing convective inertia can also be ignored in the 
first approximation, and the momentum equation becomes 

w = <^-«' 

Both the continuity (5.1) and momentum (5.5) equations are now linearized. 

In view of (5.9) the boundary condition on the shoreline (5.6) can be expressed, 
instead, as 

h^ = on S. (5.10) 
on 

Consistent with the linearized approximation, the shoreline position can be prescribed 
a priori. 

Equations (5.8) and (5.9) can be combined by the process of cross differentiation. 
First differentiate (5.8) with respect to t, 

a (DC 



then take the divergence of the product of (5.9) and h, 

v-|/i^l = -v(^wc). 



9t , 

The difference of these two equations gives 



V-(^WC). (5.11) 



For a horizontal bottom h = constant. 



1 d^C 
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where c = ^/gh — 0{\/T) is the characteristic velocity of infinitesimal wave motion. 
Equation (1.6.12) is the two-dimensional extension of the wave equation. If, furthermore, 
all conditions are uniform in the y direction, d/dy = 0, (5.12) reduces to the familiar 
form 

1 d^C 



(5.13) 



6 Sound in fluids 

The basic equations governing an inviscid and compressible fiuid are as follows. Mass 
conservation: 

|^ + V-(pu) = (6.1) 

Momentum conservation: 

p(^^ + u-Vu^ =-Vp (6.2) 
We must add an equation of state 

p^p{p,S) (6.3) 

where S denotes the entropy. When no temperature gradient is imposed externally and 
the gradient of the fiow is not too large, one can ignore thermal diffusion. The fiuid 
motion is then adiabatic; entropy is constant. As a result p — p{p, So) depends only on 
the density. Eq. (6.1 ) can be written as 

J^+„.V„]=-m V, (6.4) 



dt J \dp J 



s 



We shall denote 



so that 



C-J{^) (6.5) 



dp J 



s 



p (^^ + u • Vu ) = -C'Vp (6.6) 



It is easy to check that C has the dimension of velocity. 
From thermodynamics we also have 
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where T is the temperature and 7 = Cp/c„ = ratio of specific heats. 
For a perfect gas the equation of state is 



p = pRT (6.8) 



where R is the gas constant. Hence for a perfect gas 

dp 



= -iRT (6.9) 

^PJ s 



Liquids are much less compressible. One usually writes the equation of state as 



dp^i'A iP^'Z) ^ (6.10) 



Denoting 



1 fdp\ 



as the coefficient of thermal expansion and 

as the coefficient of isothermal compressibility. Usually j3 is small and k much smaller. 
Under isothermal conditions it is k that counts. 

The simplest limit is the case where the background density po and pressure do are 
uniform, the fluid is at rest and the dynamic perturbations are infinitesimally small. We 
can write 

P^Po+p': p^Po + p' (6.13) 
with p' <^ Po and p' <^ Po, and linearize the equations to 

^ + PoV-u = (6.14) 

and 

| = --Vp' (6.15) 

dt Po 

Taking the curl of the second, we get 

d 

—V X u = (6.16) 
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thus the velocity field is irrotational if it is so initially. We can introduce a potential 

by 

u = V(/) (6.17) 



It follows from the momentum equation 



P = -Po 



Using these we get the wave equation, 
where 



90 
di 



dpo 
dpo/ s 



(6.18) 



(6.19) 



(6.20) 



7 Flexural waves in a beam on an elastic foundation 



Referring to Figure (7), we first describe the momentum conservation of a thin beam. 
Assume the beam to have a horizontal axis when it is not loaded and to have uniform 

z\ 

\U 




Figure 6: Defiection of a beam 

material properties. Let V{x) denote the upward deflection of the beam axis. If the 
thickness is small compared to the length and the deflection small compared to the 
thickness, a plane cross section remains approximately plane after deformation. Hence, 
the longitudinal displacement U at section x and height z above the axis is proportional 
to z and to the tilt angle dV/dx 

dV 

T T r>u 
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figure 



p dx 




kVdy 



S+^dx jM+^dx 
dx J 5x 



Figure 7: Forces and moments on a beam segment from a; to a; + dx. 



Consequently, tlie strain is 



and the longitudinal stress is 



E 



dU 
dx 

dU 



dx'^ 



-Ez 



(7.1) 



ax ax^ 

where E is Young's modulus. The total moment about the mid-section ^ = due to 
the stress distribution across the section is 

52t 



-ft/a j-h/1 

M^ - I a^zdz = E— / z^dz = EI- 



-L 



h/1 



dx^ 



h/2 



dx'^ ' 



(7.2) 



where 



h/2 



z dz 



h/2 



is the moment of inertia of the cross section with respect to its mid-section z — 0. Con- 
sider a length element of the beam from x to x + dx, The balance of angular momentum 
about the center of the element requires that 

M+—dx-M+^S+ -dx^ Y + ^ Y = ^^^^a^' 

where p is the mass and pJ is the rotatory moment of inertia per unit length of the 
beam. For a beam with rectangular cross section, J = /V2, where h is the height of 
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the beam. Thus, 



On the other hand, balance of vertical forces requires that 

dS d'^V 
S + — dx — S = pdx + kV dx + p 



dx ^ ' ' ^9^2 ' 

where k is the elastic constant of the lateral support, and p(a;, t) is the distributed load. 
Making use of (7.3), we get 

which is a fourth-order partial differential equation, derived first by Lord Rayleigh. If the 
wavelength L is much greater than the beam height h/h <^ 1, then the term representing 
rotatory inertia is negligible and (7.4) may be simplified to 

Now the boundary conditions. For a beam of finite length, each end can be free, clamped, 
or supported on a hinge. At a free end, there is neither torque nor shear 

At a clamped end, the deflection and slope must vanish 

dV 

V = 0, - = 0. (7.7) 
At a hinged end both the deflection and the torque are zero 

V = 0, g^ = 0. (7.8) 

For an inflnitely long beam, the boundary conditions at inflnity depend on the loading. 
For a transient loading with flnite duration, V should vanish at inflnities. For time- 
harmonic loadings, the disturbance should at most be outgoing waves. 
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8 Homework No.l 

1. A membrane is kept taut over an area S. The membrane has a constant density p per 
unit area and is under uniform tension in all directions. Derive the governing equation 
for the lateral displacment u{x, y, t) of the membrane vibrating under distributed loading 
of p{x, y, t) per unit area. 

2. Consider the longitudinal vibration of a cylindrical rod with one end at a; = 
fixed and the other end at x = L attached to a mass M. Before t = the rod is 
compressed by the length eL with e <^ 1. At t = the compression is released. State 
the goverining equation and all boundary and initial conditions. 



3. Condsider the torsional vibration of a cylindrical rod of circular cross section of 
radius a. Let 9{x, t) — angular displacement of the cross section at x, da — area element 
in the cross section and located a the distance r from the axis, see figure 8. Let r be 
the shear stress, G the shear modulus of elasticity , and the angular displacment of a 
line orginally parallel to the axis. Show that 

Invoke Hooke's law r = G4> and show that the total torqure aplied to the cross section 
at X is 



where 

J 



J J r^da (8.3) 

is the polar moment of inertia of the cross section. Let / be the moment of inertia per 
unit length of the rod. Show that 

d^e a J d^e 



dt^ I dx^ 



(8.4) 



4. (a). Derive equations (6.1) and (6.2) for a compressible fluid. Use differential 
argument. For small-amplitude oscillations these equations can be linearized for exam- 
ining the propagation of sound. Sound waves in air and water are different because of 
the difference in equations of state: p = p{p). 



1.7. FLEXURAL WAVES IN A BEAM 
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Figure 8: Torsion of a circular cylinder 



(b). During an earthquake, water in a reservoir exerts hydrodynamic pressure on a 
dam that may fail. Formulate the dam-reservoir interaction problem under the following 
idealizations. The reservoir is infinitely long and has a uniform rectangular cross section. 

Water is present only on one side of the dam [x > 0) and has the constant depth h. 
Before t = 0, all is calm. After t = the dam is forced to vibrate horizontally so that 



The free surface is exposed to constant atmospheric pressure. The reservoir bottom is 
rigid and does not vibrate vertically (!!!). Neglect gravity but consider compressibility 
of water because of the high frequency (~ O(lOO)Hz). Express all governing equa- 
tions including the boundary conditons in terms of the velocity potential defined by 
(m, V, w) = V<p. 




(8.5) 



